Hund's coupling and the metal-insulator transition in the two-band Hubbard model 
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The Mott-Hubbard metal-insulator transition is investigated in a two-band Hubbard model within 
dynamical mean-field theory. To this end, we use a suitable extension of Wilson's numerical renor- 
malization group for the solution of the effective two-band single-impurity Anderson model. This 
method is non-perturbative and, in particular, allows to take into account the full exchange part of 
the Hund's rule coupling between the two orbitals. We discuss in detail the influence of the various 
Coulomb interactions on thermodynamic and dynamic properties, for both the impurity and the 
lattice model. The exchange part of the Hund's rule coupling turns out to play an important role 
for the physics of the two-band Hubbard model and for the nature of the Mott-transition. 



I. INTRODUCTION 

Many materials with open d- or /-shells show metal- 
insulator transitions which are commonly classified to 
be of Mott-Hubbard type due to strong electron-electron 
correlations^. The conventional model to study this type 
of transition is the one-band Hubbard model^, which in 
standard notation reads 
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Major progress in understanding the physics of the 
Mott-Hubbard metal-insulator transition (MHMIT) of 
the model has been achieved in the last decade 
through the development of the dynamical mean-field 
theory (DMFT)2d£. At T = the MHMIT occurs at a 
critical value of the Coulomb parameter U c « 1.5W^£ii, 
where W denotes the bandwidth of the density of states 
at U = 0. Interestingly, the transition is of first ordeal 
for T > with a second order end point at a T c ss 0.0171^ 
and U c « 1.2W. 

Such a second order end point is also seen in the phase 
diagram of typical Mott-Hubbard systems like V2O3 9 . 
Therefore, the one-band Hubbard model has been fre- 
quently u sed a s a microscopic model for these materials 
(see Refs. I5ll0l and, in particular, Ref. ^ which focusses 
on the critical regime close to the second order end point). 
On the other hand, the justification to base the micro- 
scopic description on a one-band model (see Ref. fl^h has 
been questioned recently^. 

For a proper description of materials such as transition 
metal oxides, the orbital structure of the relevant elec- 
tronic degrees of freedom has to be taken into account. 
This can lead to a fairly complicated form of the underly- 
ing tight-binding bandstructure (the kinetic energy term 
acquires a matrix structure). Furthermore, additional 
local Coulomb matrix elements arise which describe the 
interactions between electrons in different orbitals. 

The simplest possible extension of the model (Q to 



the case of orbital degeneracy is the two-band Hubbard 
model. It is a relevant model whenever the electronic 
degrees of freedom close to the Fermi level are two-fold 
degenerate, as for the e" g states in materials like LaMn.03 
or KCuF^i. Here we investigate the two-band Hubbard 
model in the following form: 
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with the orbital index m = 1,2. The Coulomb parame- 
ters V > and J > describe the inter-orbital Coulomb 
interaction and Hund's exchange coupling, respectively. 
The last term in (J2J is necessary to ensure rotational 
invariance of the interaction. By virtue of this rota- 
tional invariance the Coulomb parameters are related by 
U' = U — 2 J. Generally, the hierarchy of interactions is 
U > U' > J. 

A multi-band Hubbard model as in eq. @ displays a 
Mott transition at all integer fillings (not only at half fill- 
ing as in the single-band case). The additional Coulomb 
interactions also modify the value of U c and, possibly, the 
character of the transition, as has already been investi- 
gated within the DMFT frameworli 14 i 15 i 16 i 17 i 18 . 

Various theoretical and numerical techniques have 
been employed to investigate multi-band Hubbard 
models within DMFT. The Quantum Monte Carlo 
methodi^i^, which is very successful in treating the 
multi-band Hubbard models in present applications of 
the LDA+DMFT approach^, cannot, however, handle 
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the rotationally invariant form of the interaction in (0) 
due to the sign problem (fo r recen t progress in reducing 
this sign problem, see Ref. I2fl2lh . On the other hand, 
exact diagonalizationi&, linearized DMFT— and exact 
treatments in the limit of infinite orbital degeneracyi£ 
are not able to reliably calculate dynamical properties. 

In this paper, we use Wilson's numerical renormal- 
ization group (NRG)2£ to solve the effective two-orbital 
quantum impurity problem which appears in the DMFT 
for the two-band Hubbard model. In Sec. II we start with 
some technical issues related to the NRG in the two-band 
case. Section III shows thermodynamic and dynamic 
quantities for the two-orbital single-impurity Anderson 
model, with the focus on the role of the Hund's coupling. 
In Sec. IV we discuss the two-band Hubbard model; here 
we concentrate on the simplest case, i.e. degenerate or- 
bitals and intra-orbital hopping only: t™j m — tij8 mm i. 
The model eq. (0 is studied on a Bethe lattice (mainly to 
compare our results with those from other approaches); 
the generalization to other lattices (other densities of 
states) is straightforward. The paper is summarized in 
Sec. V. 



II. NRG FOR MULTI-ORBITAL MODELS 

The use of the NRG to solve the effective quantum im- 
purity model appearing in the DMFT self-consistency has 
been extensively discussed in the literatur o 7 i 8 i 23 . Here 
we focus on the additional problems arising in the two- 
orbital impurity Anderson model, which constitutes the 
effective local model arising in the DMFT for the two- 
orbital Hubbard model. 

The two-orbital Anderson impurity Hamiltonian (in 
standard notation) is given by 

H = "S~^ er ct c- + e,]dl y ,„d,„„ (3) 
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Within the NRG approach, the quantum impurity 
problem is mapped onto a semi-infinite chain form22 with 
the impurity at the first site of the chain and the con- 
duction band written in a one-dimensional tight-binding 
form. The model in the semi-infinite chain form is then 



solved by iterative diagonalization. Therefore, the ma- 
jor obstacle in applying the NRG to multi-band models 
is the dramatic increase of the Hilbert space with each 
NRG step (in each step, one additional site of the semi- 
infinite chain is included). 

A possible solution to this problem is to use very large 
values of the discretization parameter A so that the num- 
ber of states can be reduced significantly. Averaging over 
many discretizations (the so-called "Z-trick") has been 
shown to give reliable results for thermodynamic quan- 
tities (see Ref. 1^), even for large values of A. However, 
this approach becomes at least cumbersome for the cal- 
culation of dynamic quantities which are required for the 
DMFT self-consistency. 

Here we adopt two different strategies which allow to 
use a small value of the NRG discretization parameter A 
and to keep enough states in each step so that dynamic 
quantities can be calculated reliably. The first one is to 
explicitely include the orbital quantum number in the it- 
erative construction of the basis states. This additional 
quantum number significantly reduces the typical matrix 
size, so that 0(7000) states can be kept in each NRG 
iteration with reasonable computation time and memory 
consumption on modern SMP high-performance comput- 
ers like e.g. the IBM Regatta. The price to pay is that 
one has to omit the last term in the Coulomb interaction, 
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because it explicitely breaks the orbital symmetry. How- 
ever, it turns out that this term does not influence 
the thermodynamics, and dynamic quantities arc only 
slightly affected via the multiplett structure of the Hub- 
bard bands. 

The second strategy is an asymmetric truncation 
scheme: instead of adding both orbital degrees of free- 
dom simultaneously, the Hilbert space is truncated after 
adding each orbital individually. This also leads to a sig- 
nificant reduction of the Hilbert space in the iterative 
diagonalization. 

However, the asymmetric truncation scheme does not 
guarantee that orbital symmetry is preserved during the 
NRG iterations. In fact, a slight violation of orbital sym- 
metry is observed for very small energies, typically much 
lower than the Kondo temperature. It turns out that for 
the DMFT-calculations presented in Sec. II VI both meth- 
ods give almost identical results, at least on the scale 
shown in the figures. 

Figure^shows a comparison of the local single-particle 
density of states (DOS) for the impurity Anderson model 
cq. (EJ with twofold degeneracy, particle-hole symmetry, 
and calculated with both schemes (symmetric truncation 
with orbital quantum number and asymmetric trunca- 
tion) in the Kondo limit. Here, the NRG discretiza- 
tion parameter is A = 2.5 and 3600 states were kept 
in each NRG step. The model parameters are U = 7Ao, 
J = C//100 and U' = U — 2 J. As usual, A = ttN f V 2 
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FIG. 1: Comparison of the single-particle dynamics for a NRG 
calculation (orbitally degenerate and particle-hole symmetric 
single- impurity Anderson model) in the Kondo regime using 
the conventional truncation scheme (full line) and the asym- 
metric truncation scheme (dashed line). Note that particle- 
hole symmetry and orbital degeneracy are preserved even in 
the latter scheme. 

denotes the bare hybridization width (Np is the conduc- 
tion electron DOS at the Fermi energy). Obviously, the 
asymmetric truncation (dashed line) leads to accurate re- 
sults for both the value of the low-energy scale and the 
form of the Kondo resonance. In addition, it does not 
violate particle-hole symmetry and orbital symmetry, at 
least on the scale shown in Fig. ^ 

The aymmetric truncation scheme introduced here 
might be of advantage in cases where the orbital symme- 
try is violated from the outset so that the orbital quan- 
tum number cannot be used in the calculation to reduce 
the matrix size. 



III. HUND'S RULE COUPLING IN THE 
SINGLE-IMPURITY MODEL 

Before we turn to the application of the NRG to the 
two-orbital Hubbard model let us first discuss the effects 
of Hund's coupling for the single impurity model 10. For 
simplicity we consider a conduction band with constant 
DOS, p(e) = Np, in the interval [—D,D] and choose 
D = 1 as unit of energy. The local energy td is cho- 
sen such that the model is particle-hole symmetric, i.e. 
( n d) — 2. As we will see, there is a profound differ- 
ence between the cases with rotationally invariant and 
Ising-like exchange. The two-orbital impurity Anderson 
model has already been investigated with NRG by sev- 
eral groupa^^SiSi. Here we want to concentrate on the 
influence of Hund's coupling on low energy scales and 
the possibility of a quantum phase transition. To this 
end we present results for thermodynamic and dynamic 
properties. To improve the accuracy of thermodynamic 
properties we employed Oliveira's "Z-trick"— , which al- 
lows to use a larger discretization A (we used A = 5) and 
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FIG. 2: Thermodynamic properties for a particle-hole sym- 
metric two-orbital impurity model with varying Hund's rule 
coupling J — 0, [7/100, and (7/10; dashed lines: entropy, solid 
lines: effective squared moments. The Coulomb parameter U 
and the hybridization Ao were chosen such that the system is 
in the Kondo limit. The inter-orbital Coulomb parameter is 
fixed to U' — U — 2J by rotational invariance. For comparison 
the results of a single-orbital SIAM are included. 



reduce the number of states kept (1000 after truncation 
here). 

Let us begin with thermodynamic properties of 
the particle-hole symmetric two-orbital single impurity 
model in the Kondo limit, that is for a hybridization 
width Ao = irNpV 2 much smaller than the other bare 
energy scales^. The temperature evolution of the effec- 
tive squared moment^i Me// '■— T ■ Ximp{T) and entropy 
S(T) for J = 0, (7/100 and (7/10 is shown in Fig. 
The calculations were done with the rotationally invari- 
ant exchange coupling and U' = U — 2 J, but neglecting 
the term breaking orbital symmetry. For comparison we 
also include results for a single-orbital SIAM - marked 
to = 1 in Fig.0- with the same values of U and Ao- 

For J = 0, i.e. U = U', we observe an intermediate 
"local-moment regime" with entropy S" = ln6 and effec- 
tive moment Me// = 1/3 corresponding to the six degen- 
erate states, two of them magnetic, in the atomic limit. 
As expected from general SU (N) arguments^ this local 
moment is eventually Kondo screened with an enhanced 

Kondo scale ~ (TjP =1 ) m . As a technical sidemark let 
us point out that for J = 0, unlike finite J, the asymmet- 
ric truncation does not work properly, leading to wrong 
results for S and Me// m the local moment regime. 

In the atomic limit, any finite J > leads to a spin 
triplet S = 1 as ground state with moment (S%) = 2/3 
and entropy In 3. Apparently, this situation is realized 
for intermediate temperatures for both J = (7/100 and 
J = U /10. At low temperatures, this local triplet is again 
quenched by the conduction electrons like for an ordinary 
5 = 1/2 Kondo effect. Obviously, the two-orbital system 
has a considerably reduced low-energy scale Tk, which in 
addition decreases strongly with increasing J. At present 
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we do not have a satisfactory explanation for this obser- 
vation, but believe that it is related to the problem how 
spin- 1/2 electrons screen a true 5 = 1 object. 

A direct consequence of this substantial reduction of 
the low-energy scale for the application to the Hubbard 
model is that critical interactions for an MIT are also 
strongly reduced for finite J (see Sec. IIVI below). 
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FIG. 3: Local DOS at T = for a particle-hole symmetric 
two-orbital impurity model with Hund's rule coupling J = 0, 
(7/100 and 17/10. Other parameters are the same as in Fig. El 

The local DOS at T = for J = (7/100 and J = U/10 
is shown in Fig. The calculations were done with a 
discretization parameter A = 2.5 and 6400 states kept 
after truncation. On the scale used in the main panel 
of Fig. |21 the Kondo resonance for J = 17/10 appears to 
be a vertical line, pointing to a strongly reduced Kondo 
temperature, too. In addition, new structures on the 
scale of J appear as shoulders in the DOS. 

If one zooms into the region [— 10 -7 , 10 -7 ] around the 
Fermi energy (see inset to Fig. [3J| , only the resonance for 
J = U /10 remains visible with an energy scale well below 
10~ 7 . This again confirms the result from the thermo- 
dynamic quantities, viz that with increasing J an expo- 
nential reduction of the Kondo temperature occurs. It is 
quite obvious, that such a reduction in Tk will later man- 
ifest itself in a corresponding reduction of the critical U 
for the Mott-Hubbard metal-insulator transition within 
the DMFT. 

Let us point out that Friedel's sum rule implies as usual 
the constraint 7rAo7l(0) = 1. This constraint is fulfilled 
with high precision due to the calculation of the DOS via 
the self-energy according to Ref . 123 - 

A completely different picture is obtained for an Ising- 
like exchange interaction in model @, which is realized 
by replacing 
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In this case, the atomic ground state is doubly degenerate 
and consists of the two states where two electrons with 
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FIG. 4: Comparison of thermodynamic properties for a 
particle-hole symmetric two-orbital impurity model with 
Ising-like Hund's rule coupling J = (7/100 and two values 
of the hybridization parameter Ao. Apparently, there exists 
a critical Ao, c which separates a strong coupling fixed point 
(Ao > Ao, c ) from a local moment like behaviour (Ao < Ao, c ). 

the same spin occupy different orbitals. In contrast to 
the full exchange we find that the properties change quite 
dramatically with the ratio J/Aq, where Ao denotes the 
hybridization width. In Fig. 0] we compare calculations 
for J = (7/100 and two different values of Ao. For large 
Ao we find the expected screening and corresponding for- 
mation of a Fermi liquid at low temperatures. However, 
for small Ao, this behaviour is replaced by the formation 
of a state reminiscent of a local moment with entropy 
In 2 and effective local moment Me// ~ 1 ■ Obviously, the 
neglect of the spin-flip terms in Hund's exchange leads to 
a "critical" ratio J/Ao jC separating strong-coupling from 
local moment behavior. We interpret this feature in the 
following way. For the full exchange interaction, spin-flip 
scattering as in the conventional 5 = 1/2 case is presum- 
ably leading to a scenario similar to the standard Kondo 
effect. For Ising-like Hund coupling J, on the other hand, 
the atomic ground state consists, as already mentioned, 
of the two states where two electrons with the same spin 
occupy different orbitals. Quite apparently, these two 
states cannot be connected by low-energy processes like 
spin-flips, i.e. the mechanism leading to the Kondo effect 
is not present here. However, if the coupling to the band 
states is large enough such that the 5=1/2 Kondo tem- 
perature is larger than J, the system can screen the spins 
for each orbital individually before the coupling J locks 
the system into the states with S z — ±1, leading to the 
observed strong-coupling behavior at large Ao- 

Presently, it is not clear whether the change in the im- 
purity properties is connected to some quantum critical 
behavior like in the pseudo-gap model 30 . The clarifica- 
tion of this question is of course of some interest in its 
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own right and will be discussed in detail in a forthcoming 
publication. 
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FIG. 5: Comparison of the single-particle dynamics for Ising- 
like Hund's rule coupling for hybridization strengths above 
and below the critical Ao, c - Parameters are the same as in 
Fig.H 

The impurity DOS corresponding to the two different 
regimes is shown in Fig. [3J For Ao > Ao. c (full line) the 
typical structure is obtained. Decreasing A below Ao, c 
completely changes the structure of the DOS. Instead of 
a Kondo peak, we now find a structure with a pseudo- 
gap at the Fermi energy. It is quite evident, that this 
"criticality" in the impurity model has profound effects 
on the MIT in the DMFT calculations for the Hubbard 
model. 



IV. MIT IN THE TWO-ORBITAL HUBBARD 
MODEL 



In particular, a qualitative change from continuous for 
J = to discontinuous for any finite J has been observed 
in Ref. (see also the Gutzwiller results in Ref. 136^ . A 
significant quantitative change of U c when excluding the 
spin-dependent part from the exchange coupling in eq. |2 
has already been mentioned in Ref. but detailed re- 
sults have not been published yet. Recently, the issue of 
possible orbital- selecti ve Mott transitions has been inves- 
tigated in Refs. I37l38l 

Let us now discuss the results from the NRG for the 
particle-hole symmetric case. To allow a direct compar- 
ison with earlier results, we use as non-interacting DOS 
the semiclliptic form po{uj) = ~vl — uo 2 (Bcthe lattice) 
with the same bandwidth for both orbitals. As NRG dis- 
cretization paramter we choose A = 2.5 and keep 6400 
states after truncation. 

Except for J = 0, all the following results were ob- 
tained with the asymmetric truncation scheme intro- 
duced in Sec. II including all Coulomb interactions from 
the model (J2J. We have observed that the symmetric 
truncation scheme (taking into account the orbital quan- 
tum number and negelcting the last term in the inter- 
action) leads to identical results for spectral functions, 
apart from some weak redistribution of spectral weight 
in the Hubbard bands due to a different atomic multi- 
plett structure. However, the latter effect is only barely 
noticable due to the broadening introduced in calculat- 
ing continous spectra from the NRG. More important, 
the critical values U c for all two-band calculations within 
DMFT are not affected; possible problems due to the 
breaking of orbital symmetry by the asymmetric trun- 
cation scheme do not play a role here, except again for 
J = 0. 

We begin by comparing results for the one-band Hub- 
bard model (Fig. |SJi) and a two-orbital Hubbard 
model with J = and U = U' (Fig. EJ>). As is well- 



Theoretical investigations of multi-orbital Hubbard 
models within DMFT have already led to a bet- 
ter understanding of various issues such as the na- 
ture of the Mott-transition as a function of or- 
bital degeneracyi2*i£*iliiiLi2i2i and the structure of the 
spectral function in realistic treatments within the 
LDA+DMFT approach-^. Detailed results have been 
obtained for the dependence of the critical interaction 
strength U c on the number of orbitals M and different 
integer fillings n. Numerical DMFT-QMC calculations 
have been performed for M < 3iii^S. Remarkably, in 
the limit of large orbital degeneracy M — > oo an ana- 
lytical treatment of the DMFT becomes possible for the 
MIT 17 . For T = 0, a scaling U c = U c , 2 oc M for the 
actual transition is found while U c \ cx \f~M is obtained 
for the critical interaction where the insulating solution 
breaks downii. This is consistent with th e linear de- 
pendence for large M found in Refs. bs 33 .'? 1 and with 
the square-root dependence reported in Ref. I35L The 
inclusion of the Hund's rule exchange coupling J has 
been shown to significantly reduce the value of 
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FIG. 6: (a) Development of the DOS as function of U for (a) 
one-band model eq. Jfl and (b) two-orbital model eq. |JSJ with 
J = and U = U' . 

known from earlier DMFT calculations^^, the critical 
Coulomb parameter, U c , for the model with J = in- 
creases strongly with the orbital degeneracy M, U c ~ 
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M 17 ' 18 . This is also apparent from the results in Fig.|Bp. 
The system stays metallic up to the largest £7 shown. 
The actual MIT occurs for a value 2.5W < U c < 3W. 

The results in Figs. |H| and [7| are calculated with a 
fairly large value of A = 2.5 and broadening parame- 
ter b = 0.8 (see eq. (8) in Ref. |8j) for both the single- 
band case in Fig. and the two band case in Figs. HJj 
and [7| We therefore expect that our critical values £7 C 
for the Mott transition differ from more precise calcula- 
tions, and indeed we find them to be somewhat overes- 
timated. In the single-band case, for example, our result 
for U c ~ 1.6W is slightly (10%) larger than the well es- 
tablished value U c = 1A7W 7 . A similar overestimation is 
also present for M = 2, where U c 2.5W is reported in 
literaturei^i^iiSiiLiS. However, the qualitative features 
of the transition are not altered by this overestimation. 

The influence of Hund's coupling on the development 
of the spectra and the occurence of the MIT can be seen 
in Fig. UJi, where results for different values of U and 
a full Hund's exchange J = £7/4 are presented^. As 
has been noted befor o 16 i 18 , finite J substantially reduces 
U c - Such a behavior is also seen in Fig. EJu The critical 
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FIG. 7: Development of the DOS as function of U for the 
two-orbital model with (a) full J = £7/4 and (b) two-band 
model with Ising-like J = [7/4. In both cases U' = U — 2J 
was used. The assignment between the different line styles 
and values of U is the same as in Fig. 

Coulomb parameter is reduced from U c ~ 3W for J = 
to U c « 1.1W. 

Presently, the standard technique to solve quantum im- 
purities with orbital degeneracy is Quantum Monte-Carlo 
(QMC). However, due to the minus sign problem, one 
has to restrict the Coulomb interaction to density-density 
type only, i.e. an Ising-like Hund's exchange. The last 
term in J5J has to be neglected completely. This raises 
the following questions: What are the consequences of 
this approximation for the dynamics and in particular 
the MIT? 

To answer this question (at least partially) we per- 
formed calculations with Ising-like exchange interaction 
as defined in the previous section (see Ref. lioT) The re- 
sults are shown in Fig. [7p. At a first glance, the results 



are not very different, except that the critical U is fur- 
ther reduced to U c ~ 0.8T4 7 . On the other hand, the 
results for the impurity calculation in Figs. 0] and \5\ al- 
ready indicate that the replacement of Hund's exchange 
by an Ising-like term has more severe consequences than 
a mere quantitative change of energy scales. In the fol- 
lowing we show that this approximation indeed leads to a 
qualitative change in the physical properties of the Mott- 
Hubbard MIT. 

Let us now turn to the nature of the Mott transition. 
For a one-band model, it is now commonly accepted that 
the transition is of second order at T — with a quasipar- 
ticlc weight that vanishes smoothly as one approaches U c - 
There is, however, a substantial region below U c , where 
the insulator is metastable^i Previous work using the 
so-called linearized DMFT (L-DMFT) suggests that for 
orbitally degenerate systems with finite J this may be 
different^ 8 ". The authors of Ref. found a first order 
transition for small to intermediate J signalled by a jump 
in the quasiparticle weight at U c . The NRG results for 
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FIG. 8: Inverse effective mass as function of U/W and J = 
17/4 for the full (circles) and Ising-like exchange interaction 
(squares). The latter shows a strong jump in 1/m* at U c , 
leading to a first order transition, while the former vanishes 
continuously. The triangles represent a calculation with fixed 
J = 0.1. The lines are meant as guide to the eyes. 

the inverse effective mass (= quasiparticle weight) for the 
case J = 17/4 are shown in Fig. [S] The circles were ob- 
tained from calculations using the full interaction, while 
the squares represent calculations with Ising-like inter- 
actions. Apparently, the latter signal a strong first or- 
der transition at U c ~ 0.81F, while the former lead to a 
continuously vanishing quasiparticle weight. For a fixed 
J = 0.1 (triangles in Fig.JSJl the quasi-particle weight near 
U c also shows a jump at U c as predicted by L-DMFT 18 . 
However, the magnitude of this jump comes out much 
smaller in our calculations. 

The differences between the different calculations (J = 
0, J = 0.1 and J = £7/4) become more apparent when 
one looks at the local squared moment (Si). For £7 — > 
this quantity has the value (Sf) = 1/4 (for J = 0.1 it 
is actually slightly larger), while deep in the Mott in- 
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FIG. 9: (Sz) for the two-band Hubbard model and different 
values of J. The lines are meant as guide to the eye. Note 
that for J = 0.1 and Ising-like J = U/4 a discontinuity in 
(S%) occurs at the critical U. 



sulator it acquires the atomic value enforced by Hund's 
coupling, i.e. (Si) = 1/3 for J = 0, (Sf) = 2/3 for finite 
full J and (Si) = 1 for Ising-like J. This behavior is 
readily found in the calculated values of (Si) in Fig. [3] 
In accordance with the results presented in Ref. 0, the 
limiting value for (Si) is approached smoothly for J = 0. 
The same holds for J = U/4, consistent with the results 
for the quasi-particle weight in Fig. [SJ The slope, how- 
ever, strongly increases when one approaches U c from 
below. For constant J = 0.1, the numerical results are 
not decisive, and could be interpreted as both a small 
discontinuity at U c and a continuous approach with di- 
verging slope. On the other hand, a quite strong discon- 
tinuity at U c appears for an Ising-like Hund's coupling 
J = U /4, signalling a rather strong first order transition 
in this case. The above results are in rough agreement 
with the L-DMFT predictions^, although there one ob- 
serves a first-order transition also for smaller values of 
J. 



The appearance of an unambiguous and rather strong 
first order transition for an Ising-like exchange coupling 
shows that in this case the physics underlying the Mott- 
Hubbard transition is very different from the one for the 
rotationally invariant exchange interaction. As for the 
single impurity model, we believe that a transition be- 
tween individually screened orbitals on the metallic side 
to a local moment regime enforced by the Ising coupling 
on the insulating side occurs as soon as J becomes of the 
order of the Fermi liquid scale. Depending on the details 
of the non-interacting DOS (its value at the Fermi level 
and the band width) and the values of U and J, this can 
lead in the worst case to a serious underestimation of U c 
and possibly an incorrect description of the behavior of 
physical quantities close to the transition. 



V. SUMMARY AND CONCLUSIONS 

In this paper we presented first studies of the Mott- 
Hubbard transition in a two-orbital Hubbard model 
within the DMFT at T = using Wilson's NRG. In addi- 
tion to a standard NRG implementation using the orbital 
quantum number, we proposed an asymmetric truncation 
scheme which turns out to work rather well in both the 
two-orbital single impurity Anderson model and the two- 
band Hubbard model. 

As a first interesting result, we observed that for the 
particle-hole symmetric case a finite Hund's exchange 
J > leads to a tremendous reduction in the low-energy 
scale Tk- This is in striking contrast to the result for 
J = 0, i.e. U = U', where the behavior conventionally ex- 
pected for an SU(N) Kondo model, viz T K ~ \/t^ =1 , 

is founo-2. At present, the precise theoretical reason for 
this rather unexpected strong influence of J on Tk is not 
clear. Interestingly, it is also rather different from the 
case (rid) ~ 1 where we find a mild increase of Tk with 
increasing J. Obviously, a detailed study of the physics 
of multi-orbital quantum impurity models has to be an 
important future aim. 

A completely different behavior occurs if one replaces 
the rotationally invariant Hund exchange by an Ising-like 
one. In this case, the levels of the atomic doublet with 
S z = ±1 enforced by the Ising coupling cannot be con- 
nected by Schrieffer- Wolff type spin-flip processes. Thus, 
the Kondo effect can only occur for J < T K n=1 , while 
for larger J the systems is locked into a local moment 
enforced by the exchange coupling. Note that this inter- 
pretation also implies that in the strong-coupling phase 
the spins on each individual orbital will be screened sep- 
arately, while for the rotationally invariant case a full 
5 = 1 system must be screened. The details of the quan- 
tum phase transition between strong-coupling and local 
moment phases have not yet been analyzed. However, in 
view of a possible relevance of an Ising-anisotropy in the 
presence of crystal fields, a further investigation of this 
problem is certainly interesting. 

We also applied the NRG to the two-orbital Hubbard 
model in the framework of DMFT to investigate the 
Mott-Hubbard metal-insulator transition at T = for 
the half-filled case. The major goal was here to eludi- 
cate the influence of Hund's coupling on the MIT and to 
investigate how the restriction to an Ising-like exchange 
changes the nature of the MIT. Our results are in general 
agreement with previous onesiiii^iSiiLi^ In particular, 
for finite J quantities like the effective mass or (Si) show 
diverging slopes as U /* U c , possibly even discontinuities 
as proposed by the L-DMFTi±& 

For an Ising-like Hund's exchange coupling the situ- 
ation becomes qualitatively different. As can be antici- 
pated from the behavior of the impurity model, the MIT 
is strongly first order with clear jumps in the effective 
mass and (Si). Note that the former also implies a dis- 
continuous vanishing of the quasi-particle peak at T = 
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as one reaches U c . Also the physics underlying this tran- 
sition is quite different compared to the rotationally in- 
variant case, reflecting the lack of spin-flip scattering pro- 
cesses connecting the two states S z = ±1. Thus, the 
metallic phase with Ising-like interaction will be charac- 
terized by individual screening of the spins on the two 
orbitals rather than a Kondo screening of a total spin 
5 = 1. Note that this subtlety will most likely influence 
low-energy properties on the metallic side close to U c , 
but be less important for "high-energy" properties like 
magnetic or orbital ordering. 

It is clear, that the investigations presented here are 
merely a starting point to systematically study proper- 
ties of multi-orbital impurity models or correlated lat- 
tice models within the DMFT at T — > using Wilson's 
NRG. The major advantage of this method is obviously 
its unmatched ability to handle exponentially small en- 
ergy scales and nevertheless provide reliable information 
on dynamics and thermodynamics even on high-energy 



scales. Thus, at least for two-orbital models we are now 
in a position to systematically study their physical prop- 
erties and address questions that are of fundamental in- 
terest for a realistic description of, for example, transition 
metal oxides but require local degrees of freedom beyond 
a simple one-band Hubbard model. 
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Neglect of last term in in fact appears to be non-critical. 
We therefore concentrate of the spin-flip terms in Hund's 
coupling. 

The adiabatic effective squared moment T • \imp (T) intro- 
duced by Wilson— should not be confused with the isother- 
mal quantity (S%). Kondo screening can be seen only in the 
former; the latter goes to a finite constant value as T — » 0. 



